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La presente invention conceme une mdthode et un dispositif de neutralisation 
par injection controlde de gaz, de bouchons ou accumulations de liquide au pied 
d*une portion de conduite fortement inclinee par rapport k Thorizontale ou riser se 
5 raccordant a une conduite d'acheminement de fluides polyphasiques en circulation, 
par injection controlee de gaz, tels que des hydrocarbures. 

Etat de la technique 

Afin de rendre les gisements en eaux profondes ou les champs marginaux 
suffisanrmient rentables, les compagnies p^trolieres doivent developper de nouvelles 

10 techniques d'exploitation, les plus economiques possible. Ainsi, il s'avere plus 
avantageux de transporter directement le melange diphasique, constitu6 de liquide 
(huile et un peu d'eau) et de gaz, dans une seule conduite ou pipeline vers une plate- 
forme de traitement en eaux peu profondes ou mSme vers des installations terrestres, 
pour y etre separe. On utilise a cet effet une portion de conduite fortement inclinee 

15 par rapport a Thorizontale (souvent proche de la verticale) dite « riser » par les 
specialistes, qui se raccorde a la conduite en eaux profondes. Cependant, le gaz et le 
liquide etant transportes ensemble, des phenomenes d'instabilite dans les 
ecoulements peuvent se produire dans la zone de raccordement avec le riser qui 
entrainent de graves problemes d'exploitation. 

20 En particulier, lorsque les debits de gaz et de liquide a Tentree sont faibles, la 

phase liquide s'accumule dans les points bas du pipeline et bloque le passage du gaz. 
La pression en amont augmente et fmit par expulser le bouchon de liquide vers un 
autre point bas, voire jusque dans le separateur de phase situe a la sortie. Ces 
phenomenes d' accumulation peuvent diminuer la productivity et remplir de liquide 

25 des conduites destinies a recevoir du gaz, L'un d'eux, qui est plus connu parmi les 
specialistes sous le nom de « severe slugging » a donne lieu h, de nombreuses 6tudes 
soit experimentales au moyen de boucles d'essai soit par simulation par des logiciels 
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de simulation tels que par exemple le code de simulation TACITE qui a fait Tobjet 
notamment des brevets ou demandes de brevet suivants : US 5 550 761, 
FR 2.756.044 (US 6 028 992) et FR 2 756 045 (US 5 960 187), FR 00/08 200 et 
FR 00/09 889 du demandeur. 

5 Ce phenomene de formation de bouchons va etre decrit ci-apres dans le cas 

simple illustre aux Fig.o& Ton consid^re seulement une conduite faiblement inclinee 
et une colonne ou « riser » terminee par un s^parateur destine & s6parer la phase 
liquide de la phase gazeuse. 

Le liquide s'accumule dans le point bas de la conduite et tend a bloquer le 
10 passage du gaz. Celui-ci est comprim^ jusqu'au moment ou la pression en amont 
depasse la pression due au poids du liquide accumule. Un long bouchon de liquide 
est alors pouss6 par le gaz en expansion. Sous de telles conditions, on observe un 
phenomene altematif ou le liquide bloque la phase gazeuse, puis s'evacue sous la 
pression du gaz et flnalement s'accumule pour bloquer i nouveau le gaz. 

15 Le processus periodique est done, plus pr6cis6ment, le suivant : 

Etape I : le liquide s'accumule au pied du riser et bloque le passage du gaz. La 
pression monte ; 

Etape n : le niveau superieur du liquide ayant atteint le sommet du riser, la 
phase liquide se deverse dans le separateur ; 

20 Etape in : la poche de gaz atteint le pied du riser et penetre a Tinterieur du 

riser. Le bouchon se diverse dans le separateur avec une vitesse beaucoup plus 
importante. la poche de gaz « explose »dans le riser ; 

Etape IV : lorsque la poche de gaz atteint le sommet du riser, la pression au 
pied de la colonne est minimale. Le liquide retombe le long de la parol du riser. E 
25 s'accumule k nouveau au pied du riser et un nouveau cycle commence. 



ler depot 

3 



Une technique connue des specialistes sous le nom de « gas lift » permet de 
remedier a ce phenomene. EUe consiste essentiellement h injecter en pennanente du 
gaz k la base du riser pour empScher raccumulation de liquide en bas de la colonne. 
A defaut de bien contrSler le ph6nomene, on est conduit k injecter le plus souvent de 
5 grandes quantites de gaz, ce qui rend necessaire Tutilisation de moyens de 
compression importants. De plus, Tinjection de grandes quantites de gaz modifie 
sensiblement le rapport volumdtrique du jgaz au liquide (GOR), ce qui complique les 
operations de separation des phases en haut du « riser », 

La methode selon rinvention 

10 La methode selon I'invention a pour objet, par le biais d'une mod61isation des 

ph6nom6nes d'instabilite 6voqu6s ci-dessus, d'exercer un controle dynamique 
efficace sur la pression du gaz k injecter dans des conduites de fa^on k reduire au 
mieux ces phdnom^nes. 

La methode selon F invention permet de neutraliser par injection controlee de 
15 gaz, la formation de bouchons ou accumulations de liquide au pied d'une portion de 
conduite fortement inclin^e par rapport k Thorizontale ou riser se raccordant k une 
conduite d'acheminement de fluides polyphasiques en circulation. Ce controle est 
exerce essentiellement en injectant a la base du riser un volume de gaz sensiblement 
proportionnel k la variation en fonction du temps du debit massique de la phase 
20 gazeuse des fluides en circulation, et de preference sensiblement egal, quand cette 
variation est positive. 

Suivant un autre mode de mise en oeuvre, on realise le controle en modulant 
aussi le volume de gaz injecte d'une quantity proportionnelle a la variation du debit 
massique de la phase liquide des fluides en circulation, mesuree egalement au mgme 
25 pas de temps. 
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L'injection est realisee a un instant t quelconque, d'apres la variation en 
fonction du temps du debit massique de la phase gazeuse des fluides en circulation, 
mesuree a un pas de temps precedent. 

Le dispositif de mise en oeuvre permet de neutraliser par injection controlee 
5 de gaz, la formation de bouchons ou accumulations de liquide au pied d'une portion 
de conduite fortement incline par rapport k Thorizontale ou riser se raccordant a une 
conduite d'acheminement de fluides polyphasiques en circulation. II comporte des 
moyens dMnjection de gaz connectes h, la base du riser, des moyens de mesure du 
debit de la phase gazeuse des fluides en circulation, et un calculateur pour 
10 conmiander Tinjection par les moyens d*injection, d*un volxmie de gaz sensiblement 
proportionnel a la variation en fonction du temps du debit de la phase gazeuse des 
fluides en circulation, quand cette variation est positive. 

Le calculateur est adapts par exemple k commander Tinjection, par les 
moyens d'injection, d*un volume de gaz sensiblement egal a la variation en fonction 
15 du temps du debit de la phase gazeuse des fluides en circulation. 

Suivant un mode de realisation, le dispositif comporte dgalement des moyens 
de mesure du ddbit de la phase liquide en circulation dans la conduite, le calculateur 
6tant adapte a moduler le volume de gaz injecte d'une quantite proportionnelle a la 
variation du d^bit mesure de la phase liquide. 

20 Breve description des figures 

- la Fig.l montre schematiquement un pipeline faiblement incline se raccordant a 
un riser tres fortement inclin6 par rapport a Thorizontale ; 

la Fig.2 montre schematiquement un ecoulement stationnaire avec penetration 
continue de gaz dans un separateur gaz-liquide a Textremite superieure d'un 
25 riser ; 
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la Fig.3 montre schematiquement la formation d'un bouchon de liquide au pied 
du riser ; 

la Fig.4 montre le stade ou le bouchon de liquide atteint le s6parateur en haut du 
riser ; 

la Fig.5 montre schematiquement la poche de gaz penetrant dans le liquide 
accumule dans le riser ; 

la Fig.6 montre le stade ult6rieur oil le liquide retombe vers la base du riser ; 

- la Fig.7 montre la courbe de pression dans le pipeline au pied du riser durant le 
cycle pr6c6dent ; 

- les Fig. 8 A a 8E montrent schematiquement et respectivement, dans un riser, un 
mode d'ecoulement dit a buUes (8A), un mode d'ecoulement dit intermittent 
(8B), un mode d'ecoulement dit bouillonnant (8C), un mode d'dcoulement dit 
annulaire (8D) et un mode d'ecoulement dit annulaire k brins (8E) ; 

-les Fig. 9A a 9G montrent schematiquement et respectivement, dans un pipeline, 
un mode d'ecoulement dit stratifie (9 A), un mode d'ecoulement dit stratifie a 
vagues (9B), im mode d'ecoulement dit annulaire k gouttelettes (9C), un mode 
d'ecoulement dit disperse k bulles (9D), un mode d'ecoulement dit intermittent 
(9E), un mode d'ecoulement dit a petites poches (9F), et, un mode d'ecoulement 
dit k bulles allongees (9G) ; 

la Fig. 10 montre un exemple d' evolution de la pression au pied d'un riser de 14 
m de haut, sans aucun contrdle par injection de gaz ; 

les Fig. 11 A, IIB montrent respectivement les variations du debit de liquide et de 
gaz en sortie de riser, egalement sans controle ; 

les Fig. 12 et 13 A, 13B correspondent aux Fig. 10 et 11 A, IIB pour un riser de 
250m de haut ; 
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- - la Fig 14 montre revolution de la fraction volumique de liquide en sortie de 
riser sans controle ; 

- la Fig. 15, montre respectivement un exemple de variation de la pression au pied 
du riser de 14m de hauteur avant Tinjection de gaz et la stabilisation obtenue par 
5 r injection controlee en accord avec la methode (a partir de 500s) ; 

les Fig.l6A et 16B montrent respectivement dans le meme riser les courbes de 
variation en fonction du temps du debit liquide et du debit gazeux, avant 
d'exercer le controle et apres (a partir de 500s egalement) ; 

- la Fig 17 montrent dans un riser de 250 m de haut, les courbes de variation en 
10 fonction du temps de la pression au pied du riser, sans controle (traits ou points 

en gras) et avec controle (traits pointiUes) ; 

- les Fig 18 A, 18B montrent respectivement, dans le mSme riser, les courbes de 
variation en fonction du temps de la pression de liquide et de gaz, sans controle 
(traits ou points en gras) et avec controle (traits pointiUes) ; 

15 - la Fig. 19 montre les Evolutions comparees de la pression au pied d'un riser de 14 
m avec contr61e de la pression d*injection en fonction de et (traits pleins) 

et seulement en fonction de (traits pointiUes) ; 

les Fig.20A, 20B montrent respectivement, dans le meme riser, les courbes de 
variation en fonction du temps du debit de liquide et de gaz avec controle de la 
20 pression d' injection en fonction de q^ et q^ (traits pointiUes) et seulement en 

fonction de q^ (traits pleins) ; 

la Fig. 21 montrent dans un riser de 250 m de haut, les courbes de variation en 
fonction du temps de la pression au pied du riser avec controle de la pression 
d' injection en fonction de q^ et q^ (traits pointilles) et seulement en fonction de 
25 q^ (traits pleins) ; 
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- les Fig. 22A, 22B montrent respectivement, dans le meme riser, les courbes de 
variation en fonction du temps de la pression de liquide et de gaz, avec controle 
de la pression d'injection en fonction de et q^^ (traits pointill6s) et seulement 

en fonction de q^ (traits pleins) ; 

5 - les Fig.23A, 23B montrent respectivement les variations respectives en fonction 
du temps du debit de liquide et du d6bit de gaz en sortie de riser ; et 

- la Fig.24 montre schematiquement un mode de realisation du dispositif de mise 
en oeuvre permettant de neutraliser la formation de bouchons. 

DESCRIPTION DETAILLEE 

10 I-l Mise en place de la modelisation 

La modelisation des phdnomdnes d'^coulement dans le systeme form6 de la 
conduite et du riser de la figure 1 est obtenue en faisant les hypotheses suivantes. 

On se place dans une approximation mono dimensionnelle oii Ton considdre 
les moyennes des diff6rentes grandeurs sur une section transversale (droite) du 
15 pipeline ou du riser. Comme cette approximation n*est acceptable si Ton suppose le 
rayon de courbure du pipeline infmi et son diamfetre constant, la modelisation va 
porter sur les parties de part et d' autre du coude de raccordement. 

On utilise une modelisation du type flux de derive, oti Ton a une Equation de 
conservation de la masse par phase et ou les equations de conservation des quantit6s 
20 de mouvement du liquide et du gaz sont additionn^es en une seule et unique Equation 
dite Equation de conservation de la quantite de mouvement du melange. Pour fermer 
ce systeme, on va choisir une relation traduisant une loi de glissement entre les 
phases. 

Nous faisons egalement les hypotheses classiques suivantes : Tecoulement est 
25 isotherme, les fluides sont newtoniens, le gaz est un gaz parfait, le liquide est 
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incompressible (done sa masse volumique est constante) et il n'y a aucun transfert de 
masse entre les deux phases. 

On considere que le nombre de Mach du melange est petit de sorte que les 
ondes de pression se propagent avec une vitesse infmie au lieu d'une vitesse proche 
de la vitesse du son dans le melange. Les phenomenes a hautes frequences sont 
supprimes mais les ondes de taux de vide continuent de se propager avec une vitesse 
proche de celle du melange. Cette hypothese est renforcee par le fait que nous 
etudions le systeme « en transition » vers I'etat d' obstruction ou 6tat perturbe, c'est-a- 
dire au voisinage du stationnaire. Cette hypothese se traduit dans la modelisation par 
Tabsence des termes d'inertie dans Fequation de conservation de la quantite de 
mouvement. 

1-2 Choix des grandeurs fondamentales 

On definit les quantites suivantes : 

et sont les fractions volumiques de gaz et de liquide dans les canalisations. 

®^ 9l sont les debits massiques de gaz et de liquide par unite de section. 

et sont les vitesses du gaz et du liquide. 
Pc Pt^ont les masses volumiques du gaz et du liquide. 
P designe la pression moyenne du melange. 

Toutes ces grandeurs sont des fonctions de RxR (espace et temps), nulles sur 
RxRl et continues sur RxRl. Suivant le contexte, on precisera toujours si Ton se 
situe dans le pipeline ou dans le riser. Ainsi on travaillera au niveau du pipeline, dans 
[0,L]x/J en 6crivant les variables (x,/), et au niveau du riser, dans [-H,0]xR en 
ecrivant les variables (2,r) . 
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Toutefois, si une equation est valable aussi bien dans le pipeline que dans le 
riser, on 6crira fonnellement Tequation avec les variables (xj) G [- l]u i?^. . 

1-3 Etablissement des equations intrinseques 

Toutes ces grandeurs sont liees par des relations algebriques et differentielles 
5 qui ne dependent pas de Tecoulement considere (ces ecoulements seront decrits dans 
la section suivante), c'est ce que nous appelons les equations intrinseques. Les autres 
types d'equations, principalement les lois de glissement, seront etudiees dans la 
prochaine section. 

1-3-1 Equations algebriques 

10 Tout d'abord, exprimons les relations issues directement de la definition des 

grandeurs. 

^g^PgRoVg (LI) 
Ql^PlRlVl (1.2) 

15 L'equation des gaz parfaits permet d'etablir entre la pression et la masse 

volumique de gaz, la relation suivante : 

PoM^^^ (1.4) 



— — ou Mg est la masse molaire du gaz, Tla temperature du 

melange et R la constante des gaz parfaits. a correspond a la vitesse du son dans le 
20 gaz h, 1 bar. 
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1.3.2 Equations differentielles 

La conservation de la matiere de chaque phase impose que : 

La conservation de la quantite de mouvement suivant les hypotheses 
prec^dentes, impose que : 

dp 

— = sin 0(/?^p^ + R^p^ )- (1.7) 

oil d est Finclinaison de la canalisation, g constante de pesanteur et le frottement 
parietal (frottement du flot contre la paroi). 

A ces equations, on ajoute une Equation de fermeture, sous forme d'une loi de 
glissement algebrique de la forme : ^(P, /?c , , , ) = 0 

1.4 Loi de glissement 

La loi de glissement T choisie depend du regime de Tecoulement. On peut 
considerer trois types d'ecoulement : stratifie, disperse a buUes (voire simplement a 
bulles) et intemiittent. Tous les regimes sont illustres dans les figures 8A ^ 8E et 9A 
a 9G. Pour notre etude, on considere le cas d'un ecoulement intermittent dans le riser 
et stratifie dans le pipeline. On verra par la suite, compte tenu de simplifications 
possibles qu'aucune loi de glissement ne s'avere necessaire pour le regime stratifie. 
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I. 4.1 Regime intermittent 

Ce regime provient d'une « superposition » entre un 6coulement k buUes et un 
6coulement stratifie. Lorsque le debit de gaz augmente, les bulles se regroupent et 
entrent en coalescence. De grandes bulles en forme d'obus apparaissent. EUes sont 
separees par des bouchons de liquide qui en general contiennent de petites bulles de 
gaz. 

Dans ces conditions, la loi de glissement en intermittent s'exprime : 

Vo -CoiRoVo +i?t^J-v: =0 a.io) 

d'oii la fonction '¥-^, {Vo.V^, Rl^vJ. 

De son cote, est determinee experimentalement et est de la forme : 

V„ = (O.35sin0 + O.54cos0)Vi^ 

II ETUDE DU SYSTEME PIPELINE-RISER SANS FROTTEMENT 

Dans cette partie, nous fixons = 0 dans Tequation 1.7 compte tenu que le 
riser est vertical et que Ton fait Thypothese que dans les perturbations d'ecoulement 
observees, les frottements influent peu par rapport a la gravite. La resolution des 
Equations s'en trouve a priori simplifi^e. 

II. 1 Etat stationnaire dans le riser 

Dans le cas stationnaire, les differentes grandeurs ne dependent que de 
Tespace ; elles seront ecrites surlignees pour mettre en evidence que ce sont des 
grandeurs stationnaires. De plus, toute fonction constante sera desormais notee h 
Taide d'une etoile. Ainsi, par exemple, en stationnaire, les debits de liquide et de gaz 
9g sont constants. Ds sont done not6s et ql . En outre, on appellera la 

pression en haut du riser. 
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Ecouleraent intermittent 

Uequation de fermeture est done donnee par LIO. En outre, comme 
precedemment, dans Tetat stationnaire, les debits de gaz et de liquide etant constants, 
ils seront notes de la meme maniere. 

La formulation implicite de P en intermittent s'ecrit : 



(n.i5) 



avec g donnee par la relation 



gpIc^G + pX+Co^I) 



10 La pression est calculable sous forme implicite en regime stationnaire dans 

I'hypothese du regime intermittent. 

P(z) peut done Stre supposee connue au moins num6riquement. On peut 
alors exprimer toutes les autres grandeurs en fonction de P(2) et de z. Pour les 6tats 
transitoires, on supposera done les grandeurs stationnaires connues. 

15 II.2 Modele stationnaire dans le pipeline 

Nous avons a priori le choix des regimes. Cependant la pratique montre qu*il 
est plus rationnel de choisir le regime stratifie, ou a la limite intermittent, tout en 
sachant que le second est moins probable. La resolution du regime intermittent est 

analogue a celle de la section precedence en n'imposant pas {0 = — ) dans Tequation 

20 1-7. Le systeme est une fois de plus integrable. En prenant comme hypothese que le 
regime dans le pipeline est stratifie, Tequation de fermeture est trop compliquee pour 
permettre une resolution des equations. De ce fait, nous ne faisons intervenir le 
pipeline que par Teffet qu'il a sur le riser, c'est-a-dire en mesurant les grandeurs 
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d'entree du riser (debits, pression, fractions surfaciques) au lieu de les determiner par 
rapport au calcul des grandeurs de sortie du pipeline. On remplace done les calculs 
par des mesures d'entree du riser. 

II.3 Etude transitoire dans le riser dans le modeler lineaire tangent 

5 Nous allons ci-apres d6velopper les calculs autour de T^tat stationnaire et 

lindariser les equations intrinsfeques au premier ordre pour r6soudre simplement le 
systeme transitoire. 

II.3.1 Invariant de Riemann 

L'introduction de cet invariant sert h faciliter la resolution pertinente des 
10 equations en lineaire tangent II s'exprime sous la forme k = ^^^^ qui n*a pas de 
signification physique evidente. On a done la relation suivante : 



V(z,0 €] -H.OlxRl , 1^ + 1^ = 0 (n.28) 
ot dz 

La quantite k se transporte dans le riser, en regime intermittent, a la vitesse 



15 Ainsi, la propriety de k nous permet d'affinner que la propagation de k en 

regime intermittent s'exprime : 

k(z,t) = ki-H,t-Tizy), (n.30) 



avec nz) = jlfi:^dx 



20 



La quantite k a une altitude z du riser a la meme valeur qu'en bas du riser 
T(z), une unite de temps a plus t6t. 
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Munis de cet invariant de Riemann, nous allons exprimer F(-Hj), qui est 
la variation de pression en bas du riser ou debutent les perturbations, en fonction des 
debits d*entree du riser. Ceci pourra nous permettre de concevoir un controle sur cette 
pression. 

II.3.2 Loi d'entree-sortie en transitoire, regime intermittent 

Nous supposons ici que les valeurs des grandeurs stationnaires sont donnees 
(on connait ces fonctions grace a Tetude en mode stationnaire). En appliquant les 
techniques de linearisation, nous pouvons 6tablir le r6sultat suivant : 



(n.3i) 



10 ou N^, q^^ ,q^^ sont des fonctions de / issues respectivement de q^, q\ 

ou Ton a fix6 la premiere variable respectivement a z, -H et -H. 

En effet, la linearisation de I'equation 1,7 donne : 

-3^ = -5((Pg -p2)/fc+^GP'<^)) (11-32) 

Nous allons remplacer R'^ par son expression en fonction de P' et k\ Pour 
p R 

15 cela, notons que k = ^ — . En lin^arisant cette relation, nous obtenons que 

(l-Co)Ra 

k'd^C R ) o* R 

R'c- =r— ^ — - ^ ^ • En remplagant ensuite R'^ dans 1132 par cette 

Pg 

expression, on aboutit ^ Tequation differentielle suivante : 

^a.o^ «^°<y^^'-> rfco=- ^<^° -r^^'-^'^o^ -u.') (n.33) 
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La propriete d'invariance de Riemann de k se prolonge pour sous les 
formes : 

^ + ^^^ = 0 et k\zj)-^k'(-Hj-Tiz)) avec f(z) = f i^&c 
Grace k ceci nous pouvons 6crire : 



De plus, nous pouvons calculer k' en fonction des debits, et nous obtenons 



r /' 



kC, 



0 „' 

<1l 



En injectant cette relation dans n.33, cette Equation devient : 



10 



Pa 



9c 



q'^ i-H,t-T(z)) {-H,t -T{z)) 

P*L 

n.34 



On pose alors : 



giPa-p*L)0-2<^M k 
Pa+C^k qc 



Ensuite, on intdgre cette Equation en utilisant la m^thode de variation de la 
constante. On obtient alors : 
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expJ*^-5U)<ix 



r 



Cix) 



P\z,t) = P\-H,t) + 



q'^ i-H,t-T(z)) ^q\irH,t-T{z)) 

P*L ) 



(X — \ 

expj_^B(u)du 



dx 



Enfin on pose : 



et si on pose N{z,v) - l[o:r(2)](^)^(^) ^^^^ 



F 

plus haut, 



, comme on Ta vu 



Dans ces relations, 

q'c'H » designe la variation en fonction du temps du debit massique de la 
phase gazeuse dans le fluide polyphasique en circulation au pied du riser, i.e. & 
r altitude -H par rapport au sonmiet du riser ; 



q\-H^ la variation en fonction du temps, ^ la meme altitude, du debit 
massique de la phase liquide dans le fluide polyphasique. 
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III Mise en place du contrdle d'injection de gaz 

Le mode d*action choisi pour corriger les perturbations, est de controler la 
pression en bas du riser. En effet, si cette pression reste proche d'une valeur 
stationnaire, cela signifie que les bouchons ne se forment pas et done que le gaz n'est 
jamais vraiment bIoqu6, L*action portera done sur la pression P(-//,r)en bas du 
riser. 

Ill .1 Contrdle par gas lift 

On suppose que la pression de sortie est fixee et done que Pq = 0. Ainsi, la 
relation 111-31 peut s'ecrire : 



10 



(0 



m-42 



On voit imm^diatement que si I'un des deux membres du produit de 
convolution est nul ou tres petit, alors il en est de meme pour P'i-H,t) . Or, ne 



( C k 

peut pas 6tre modifi6. Reste la quantity \q'G „ ^q'L_f, 

y pL 

d6bit, que Ton notera Q(t). 



if), homogene h. un 
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Le principe du mode de contrdle selon I'invention va consister 
essentiellement, a intervalles d'action ddfinis, h. mesurer ^ un instant ti la quantit6 ci- 
dessus, pour obtenir une mesure Af(r,), puis, ^ I'instant r, +Af, ^ ajouter k 
Q{ti + Ar) , la quantity + Af) = -M (t^ ) et ainsi de suite. On a ainsi Vi > 1 : 



«(/, + zA/) = -M(/, + (t - 1)A0 

+ lA/) = tt(r, + / Af) + Af (f 1 + tA/) 
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Bien sflr, quand Ar — > 0 , cela revient a annuler le second terme du produit de 
convolution. Comme il n'est pas physiquement envisageable d'avoir une commande 
qui aspire le gaz dans un melange diphasique, on utilise une commande telle que : 

u(t^ + f AO = max(-M (t^ + (i - l)Ar), 0) . 
5 IV Validation des contrdles par injection de gaz 

Pour valider ce mode de controle, on a utilise le code TACITE precite qui 
simule les ecoulements de la fagon la plus precise et la plus realiste possible. 

Nous avons etudie deux cas de formation de bouchons, tres diffdrents, et test6 
sur eux notre mode de controle. 

10 IV. 1 Remarques concernant les simulations 

Le logiciel TACITE utilise une mdthode de type Volumes Finis pour simuler 
les ecoulements dans les conduites et le pipeline est discretise par exemple selon la 
methode de maillage decrite dans la demande de brevet FR EN 00/08200 du 
demandeur. Pour simuler Finjection de gaz en pied de riser, nous modifions le flux 
15 entre les deux mailles situdes juste avant et juste apres le coude, numerotees 
respectivement n - 1 et Nos simulations ont conwne etat initial le regime 

C k 

stationnaire. Le coefficient — ^ est done mesur6 une fois pour toute a I'instant 

Pl 

initial. Nous identifions les debits stationnaires aux debits d'entree de pipeline, et nos 
debits en entree de riser aux debits de la maille n - 1. Pour Taffichage des resultats, 
20 nous approximons le debit de gaz a injecter par la difference de debit entre les 
mailles n etn - L 
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IV.2 Simulations des differents cas non contrdles 

Les simulations effectuees en utilisant le code de simulation TACITE, se 
basent sur deux geometries oh L represente la longueur de pipeline, H, la hauteur du 
riser et 0, leur diametre conimun. 

5 Cas 1 : L= 60 m, H = 14 m, 0 = 5 cm, 

Cas 2 : L 1750 m, H = 250 m, 0 = 25 cm 

Cas 1 : Riser de 14 m 

Voici le cas que nous avons pris pour reference : 

Po = 1 bar, Ql = 2.10'^ kg/s et Qg = 2.10'"^ kg/s. Ce cas presente des 
10 oscillations marquees, avec une periode d'oscillation autour de la mdnute. Le regime 
d'^coulement correspond h nos hypotheses : stratifie dans le pipeline et intermittent 
dans le riser. 

Les oscillations du systeme (Fig. 11 A, UB) ne sont pas assez importances 
pour annuler compldtement le debit de gaz en pied de riser ou pour que la surface 

15 libre du liquide descende en dessous du niveau de Tembouchure du riser. Nous ne 
sonunes done pas dans un cas de formation de bouchons tres importants. Cette 
propri6t6 permet de rester proches de notre hypothese de linearis^ tangent, en 6tant 
toujours dans un regime d*6coulement peu 6loign6 du stationnaire. Nous restons done 
pour tout r, dans la situation quasi-stationnaire de transition vers la formation de 

20 bouchons. 

On distingue bien deux phases principales dans le cycle : 

1. Une phase d'accumulation du liquide dans le riser, le debit de liquide 
est nul en sortie et la pression augmente ; 



2. Une phase de production du bouchon de liquide ou la pression 
25 diminue. 
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Cas 2 : Riser de 250 m 

Les caracteristiques principales de ce cas sont les suivantes : 
Po = 10 bar, Qt = 4 kg/s et Qg = 0,5 kg/s. 

On est plus proche ici des conditions d'exploitation reelles d'un pipeline. En 
5 Tabsence de tout controle, le systeme entre dans une phase de formation d'un 
bouchon de liquide. On constate que le d6bit de liquide en sortie passe par des phases 
d'expulsion deja violentes, avant meme la phase d'expulsion du bouchon et de chute 
de la pression. A la fin de la simulation, la pression a atteint un maximum, et le 
systeme se prepare a entrer en phase d'expulsion du bouchon de liquide, Le systeme 
10 n'est done dans la phase de transition vers la formation de bouchons que dans les 
premiers instants, car la simulation prend pour point de depart le regime stationnaire. 

Sur la Fig. 14, on peut observer des cycles tres rapides dans revolution de la 
fraction de liquide. C'est un signe de la grande instabilite de ce cas. 

rV.3 Controie des deux cas 

15 Le premier contrSle que nous allons examiner est le controie theorique trouve 

dans le chapitre precedent. A la suite de nos observations, nous allons voir que Ton 
peut imaginer un autre controie, que nous testerons egalement. 

IV.3.1 Essai du controie theorique 

Nous introduisons ici un controie par le gaz lift, en injectant en pied de riser 
20 un debit de gaz de la forme : 



P L 



1er depot 



21 



Cas 1 : Riser de 14 m 

Nous commen^ons a controler un peu avant t = 500 s. Las graphiques des 
Fig.IV-16, Fig.IV-17 montrent en paralldle revolution du systdme libre (points 6pais) 
et celle du systdme contrSle (pointillds fins). Ce controle parvient a maintenir la 
5 pression et les debits de sortie dans un voisinage de leur valeur stationnaire. La 
masse de gaz injectee est petite par rapport k ce qui provient du pipeline ; 
I'augmentation du debit de sortie due a I'injection de gaz est de moins de 5%. 

Cas 2 : Riser de 250 m 

Pour cette simulation, nous commengons le controle a partir de r = 0 s. Les 
10 graphiques des Fig, 15 et Fig.l6A,16B montrent en parallSle revolution du systeme 
libre (points epais) et celle du syst&me contrdle (pointilles fins). Les debits de gaz en 
presence du contrSle sont beaucoup plus reguliers que dans le cas precedent. 

Le contrdle par injection de gaz en utilisant la formule thdorique fonctionne 
done dans les deux cas, bien qu'ils soient tres diff^rents. 

15 IV3.2 Controle simplifie 

Comme on observe qu*il y a une forte correlation entre les evolutions des 
debits de gaz et de liquide, et que Ton a presque toujours |^*^ ^ | > 
possible aussi de controler le systdme avec la commande simplifiee : 

M(0 = max(0,(9'^^^ )(r-Ar)) 

20 Cas 1 : Riser de 14 m 

Dans les graphiques des Fig.l7, Fig.lSA, 18B, le controle avec qc et qt est 
represente en trait epais, tandis que le controle portant uniquement sur qo est en 
pointilles fins. Le controle commence ici des t = Os. 



C^k 



P L 



y il est 
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Le controle portant uniquement sur qo donne des resultats pratiquement 
identiques k celui deja teste, voire 16gerement meilleurs en ce qui conceme les 
oscillations de la pression et la rapidite du controle. 

Cas 2 : Riser de 250 m 

5 Pour ce cas nettement plus instable, nous ne rencontrons pas (Fig.22A, 22B) 

d'oscillations rapides du debit de gaz avec le controle simplifie. Le systeme controle 
avec evolue avec une legere avance par rapport au meme cas controle avec 
Texpression complete. 

On constate done que le controle simplifie sans le coefficient en , permet 
10 egalement de controler le systeme mais le controle theorique est meilleur cependant 
parce quMl n'entraine pas, contrairement au contrSle simplifie, de probleme a la sortie 
du riser, et il est plus 6conomique en gaz d*injection et le systeme est aussi bien 
controle. 

rV3.3 Robustesse du contr51e 

15 Nous testons la robustesse de notre controle theorique {q^ et q^^) par rapport 

au temps de reaction des capteurs et des actionneurs. Pour cela, nous comparons le 
contrdle obtenu plus haut en ajustant le d6bit de gaz a chaque pas de temps (pointilles 
epais) avec le contrdle obtenu en ajustant ces parametres avec une frequence plus 
faible (pointilles fins). 

20 Cas 1 ; Riser de 14 m 

Ici, on reajuste le contrdle toutes les 3 secondes, ce qui correspond a 20 pas 
de temps de TACITE. Le controle conmience ici des f = 0 .y . 

On voit sur les Fig.23A, 23B, que les debits sont un peu plus irreguliers, mais 
le systeme continue a etre controle. 
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Cas 2 : Riser de 250 m 

Le systeme n'est pas stabilise avec une constante de temps de deux secondes. 
Avec une constante de temps de 1 seconde, le systeme est stabilise, mais on reste 
proche de Tordre de grandeur du pas de calcul (environ 0,3 s pour ce cas). II n'y a en 
5 outre pratiquement pas de difference avec le cas ou le controle est ajuste a chaque 
pas de temps. 

Dans tous les cas on a un bon controle meme si les debits oscillent tres 
faiblement autour d'une valeur moyenne. 

Le dispositif de mise en oeuvre comporte (Fig.24) des moyens 1 d'injection 
10 de gaz connectes a la base du riser, des moyens 2 de mesure du debit de la phase 
gazeuse des fluides en circulation, et un calculateur 3 pour commander T injection, 
par les moyens dMnjection 1, d'un volume de gaz proportionnel et de preference 
sensiblement voisin de la variation en fonction du temps du debit de la phase gazeuse 
des fluides en circulation, quand cette variation est positive. 
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REVENDICATIONS 



l)Methode pour neutraliser par injection controlee de gaz, la formation de 
bouchons ou accumulations de liquide au pied d*une portion de conduite fortement 
5 incline par rapport h Thorizontale ou riser se raccordant a une conduite 
d'acheminement de fluides polyphasiques en circulation, caracterisee en ce que Ton 
injecte a la base du riser un volume de gaz sensiblement proportionnel a la variation 
en fonction du temps du debit de la phase gazeuse des fluides en circulation, quand 
cette variation est positive. 

10 2)Methode pour neutraliser par injection controlee de gaz, la formation de 

bouchons ou accumulations de liquide au pied d'une portion de conduite 
sensiblement verticale ou riser se raccordant a une conduite d'acheminement de 
fluides polyphasiques en circulation, caracterisee en ce que Ton injecte a la base du 
riser un volume de gaz sensiblement €gal a la variation en fonction du temps du debit 

15 de la phase gazeuse des fluides en circulation, quand cette variation est positive. 

3) Methode selon la revendication 1 ou 2, caracterisee en ce que Ton module 
le dit volume de gaz injecte d'une quantite proportionnelle a la variation du debit de 
la phase liquide des fluides en circulation. 

4) Methode selon Tune des revendications precedentes, caracterisee en ce que 
20 Ton injecte & un instant t, la variation en fonction du temps du debit de la phase 

gazeuse des fluides en circulation, mesuree k un pas de temps precedent. 

5) Dispositif pour neutraliser par injection contr616e de gaz, la formation de 
bouchons ou accumulations de liquide au pied d'une portion de conduite fortement 
incline par rapport a Thorizontale ou riser se raccordant a une conduite 

25 d'acheminement de fluides polyphasiques en circulation, caracterisee en ce qu'il 
comporte des moyens d' injection de gaz connectes a la base du riser, des moyens de 
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mesure du debit de la phase gazeuse des fluides en circulation, et un calculateur pour 
commander Tinjection par les moyens d'injection d'un volume de gaz sensiblement 
proportionnel h, la variation en fonction du temps du debit de la phase gazeuse des 
fluides en circulation, quand cette variation est positive. 

5 6)Dispositif selon la revendication 5, caracterise en ce que le calculateur est 

adapte a commander Tinjection par les moyens d*injection d'un volume de gaz 
sensiblement egal & la variation en fonction du temps du debit de la phase gazeuse 
des fluides en circulation. 

7)Dispositif selon la revendication 5 ou 6, caracterise en ce qu'il comporte 
10 des moyens de mesure du d6bit de la phase liquide en circulation dans la conduite, le 
calculateur etant adapte a moduler le volume de gaz injecte d'une quantite 
proportionnelle ^ la variation du d6bit mesur^ de la phase liquide. 
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